Lattice Dynamics of the Heisenberg chain coupled to finite frequency bond phonons 
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The phonon dynamics in a one dimensional Heisenberg spin chain coupled to finite- frequency bond 
phonons is studied. We present the first detailed phonon spectra for these systems using Quantum 
Monte Carlo. The quantum phase transition is dominated by a central peak, yet the renormalisation 
of the main phonon branch depends strongly on the bare phonon frequency uo. The main branch 
remains largely unaffected at cuo^J, but it softens completely when ooo is low enough. This is 
an unusual scenario for a structural phase transition and was observable only on sufficiently large 
systems. Approaching the dimerized phase from finite temperature, the lattice dynamics mirrors 
the behavior of a three dimensional system. For the efficient measurement of Greens functions, we 
introduce a mapping from the stochastic series expansion to continuous imaginary time. 
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Due to the discovery of the inorganic spin-Peierls 
compound CuGeOs, one dimensional spin systems cou- 
pled to phonons have been studied intensively with 
analytical and numerical methods like density-matrix- 
renormalisation [l[ , exact diagonalisation [2j] , linked clus- 
ter expansion [3], continuous unitary transformation [4|, 
renormalisation group [5], and Quantum Monte Carlo 
(QMC) p. Within a spin-Peierls transition, magnetic 
interactions cause a structural phase transition. It can 
be classified into two categories: the displacive transi- 
tion exhibits a soft phonon mode whereas the order- 
disorder transition is characterized by an additional cen- 
tral peak without any phonon softening 0-The Peierls- 
active phonon mode therefore possesses either one or two 
time scales. This, however, is a more phenomenological 
classification and even paradigm examples of displacive 
phase transition, like SrTiOs or KMnFs, show an addi- 
tional central peak above the critical temperature apart 
from the soft phonon mode Within the spin-Peierls 
transition there exist prominent candidates for either sce- 
nario: Copper bisdithiolene (TTF-CuBDT) has a pre- 
cursive soft phonon mode at high temperature [9] and 
is therefore believed to belong to the soft phonon sce- 
nario. For CuGeOs, however, although the central peak 
was not found yet directly, the Peierls active phonon 
mode hardens [10] which is in agreement with the central 
peak scenario of the RPA [ll|, ll2[ . Gros and Werner [12[ 
showed that RPA predicts softening only for bare phonon 
frequencies smaller than some multiple of the transition 
temperature. 

In this paper we focus on the phonon dynamics of 
the spin- 1/2 anti- ferromagnetic Heisenberg chain inter- 
acting with finite frequency bond phonons. This system 
relates to real materials as mentioned above, in which 
the critical fluctuations are quasi one dimensional in a 
large temperature range above the critical temperature 
[131 ]. Using sum rules, Sandvik and Campbell [6[ have 
indirectly shown that the q = tt phonon exhibits a cen- 
tral peak even below the phase transition, for arbitrarily 



small spin-phonon coupling at T = 0. Therefore they sug- 
gested that the zero temperature transition is of the cen- 
tral peak type. We verify this conjecture by calculating 
the whole phonon spectrum. Contrary to the expecta- 
tions of ref. [6j], we also find softening, and a qualitative 
dependence of the phonon spectra on the bare phonon 
frequency. To our knowledge these are the first spectra 
for spin-Peierls transitions from unbiased QMC. We also 
discuss the lattice dynamics at finite temperature. 
The model under consideration is 



H- 



N 
i=l 



N 



(1) 



where xi = -j= (ai + a^j is the spatial amplitude of the 

phonon which affects the lattice bond + 1), g is the 
coupling between spins and phonons, uoq is the bare fre- 
quency of the dispersion- less Einstein phonons (EP), and 
rii = a\ai is the phonon occupation number operator cor- 
responding to Xi. We use J as our unit of energy. The 
system is studied with a Quantum Monte Carlo (QMC) 
simulation based on the stochastic series expansion (SSE) 
The spins are updated globally using the loop 
algorithm 16|, 12] whereas the phonons are updated lo- 
cally. To overcome the resulting high autocorrelations at 
large coupling g, we have used simulated tempering [lit 
in g , with an automated choice of couplings and weights 

a 

We study the phonon dynamics by analyzing the struc- 
ture factor S x (q,w). It is obtained via the inversion of 



(x- q (r)x q (0)) = / d(jS x (q,cj) e 



(2) 
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The measurement of such Greens functions in the 
SSE representation has been very cumbersome [2y]. In- 
stead, we make use of a correspondence between SSE 



where x q ~^E'-=i 
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and continuous imaginary time [l4(. For a Hamiltonian 
H = Hq — V, the interaction representation of the parti- 
tion function is 



TrT e-^» / dr n ... / dr 2 / 

n=0 J JO i0 



dTiV(Ti)...y(T„), 

(3) 

with V(r) =e HoT Ve HoT . Indeed, this is the represen- 



21, 



tation of continuous time QMC [6, 
When i?o = 0, V(r) does not depend on r, and the time 
integrals just provide a factor /3 n /n!, which multiplies 
V n . When also V = H\ + H2 + . . . , this product is a sum 
of ordered operator products, which are the SSE opera- 
tor strings. Given an SSE configuration with n opera- 
tors, we can now stochastically map back to continuous 
time by randomly choosing n times between and /3, 
ordering them, and assigning them to the n operators. 
We obtain a configuration in space and continuous time. 
The Monte Carlo ensembles thus obtained in SSE repre- 
sentation and in imaginary time are equivalent. Greens 
functions can then be measured efficiently by evaluating 
the measurements on a fine time-grid and using the Fast 
Fourier Transform in space and time. Care should be 
taken with off-diagonal operators. Details will be pub- 
lished elsewhere [19|. In our simulations, this procedure 
reduced the computational effort for measuring Greens 
functions in all of momentum space by about a factor of 
57V, where N is the system size. 

We improved the inversion of ([2j) with the maximum 
entropy technique (MaxEnt) [25] by using the exactly 

known moments M 1 = and M 3 = where M n = 

00 

J duo S (g, uj) uo n . These moments roughly double the 

-co 

resolution of S(q,ui) near uj ~ uo^. The temperature in 
our simulations was chosen low enough to achieve the 
T = limit where no changes in S (q,uu) by lowering the 
temperature are observed. We used at least /3 = 27V. 

The spin-Peierls chain Eq. (pQ) exhibits a Kosterlitz- 
Thouless transition [y, [26[ at a finite coupling g c [6[ to 
a dimerized phase with a finite spin gap A5, which be- 
gins, however, exponentially small in g — g c . Citro et al. 
[26[ identified the relevance of the ratio of As to ujo: at 
large As the system is adiabatic, while it is in an anti- 
adiabatic regime at small As/ojq. This is the case for 
finite frequency Einstein phonons until far beyond the 
phase transition. Sandvik and Campbell [6] showed that 
for 0<g<g c the phonon structure factor, which is 0.5 at 
g = 0, already diverges like log TV, reflecting a long range 
1/r behavior in the phonon correlation function. 

We determined the critical coupling g c in the same way 
as ref. [6[ from the spin susceptibility. Xs/N diverges 
for g < g c (like log^TV for the pure spin model g = 0), 
becomes constant at g = g c where logarithmic corrections 
vanish, and diminishes at larger g, reflecting the finite 
spin gap. Our results for g c agree very well with recent 
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FIG. 1: Dynamic structure factor S x (k,Lj) at small uo — 0.25 
in a density plot, for spin-phonon couplings (a) g = 0.1, (b) 
g = 0.23 close to the critical coupling, and (c) g = 0.3, at 
lattice size N — 256 and (3 — 512. The error bars indicate the 
position of peaks and their uncertainties. The insets show 
S x (tt,uj) and the MaxEnt resolution (see main text). The 
shaded blue area has a weight of 0.5, equal to the total spectral 
weight at g = 0. For better visibility, the height of the low 
energy peak (central peak) at k = tt has been truncated to 
the second highest peak in the spectra for g = 0.23 and 0.3. 
Depending on a central peak as well as softening occur. 



DMRG calculations (Fig. 1 in ref.[27|). 

We will now discuss the phonon spectra. They differ 
between small and large c^o- We will first present results 
at small cjo = 0.25 (FigJTJ). Ref. [6] indirectly deduced the 
presence of a central peak from the ratio of the phonon 
susceptibility and the static phonon structure factor. The 
latter grows like log TV up to the transition, and like N 
beyond. 
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Indeed, at low coupling g = 0.1 (FigQJi) we find a cen- 
tral peak and in addition a new phonon branch start- 
ing close to uo = and q = tt which is generated by the 
spin-phonon interaction. We will call it the central peak 
branch (CPB). The finite energy gap of the CPB at q = TT 
scales to zero as 1/N. Away from q = tt the spectral 
weight of the CPB decreases rapidly, leaving the EP well 
separated from the CPB. The dispersion of the CPB (dot- 
ted line) coincides with the spin velocity, as discussed 
later. 

The insets of FigQ] provide some insight into the dis- 
tribution of spectral weights. We use a test spectrum 
with two narrow peaks at the same position and with the 
same weight as S x (ix,ui). We add noise similar to that 
observed in (x-^ (r) x^ (0)) to the correlation function 
obtained from the test spectrum. The resulting MaxEnt 
spectrum is shown as dashed lines. In the case of FigOJi 
this shows that S(tt,u) is significantly wider than the 
MaxEnt resolution. The shaded area under S(ir, uo) has 
a spectral weight of 0.5 which is the total spectral weight 
without spin-phonon interaction. The EP at q = tt there- 
fore widen but preserve their spectral weight. 

At # = 0.23 ~ g c (FigHb) the EP are strongly affected 
by the spin-phonon interaction. They soften partially. 
The EP peak becomes broader and less well defined. 
Some spectral weight even extends to above c^ = 0.3. The 
bare phonon frequency is not reflected in the spin spectra. 
This can be seen in Fig|2] which shows the dimer-dimer 
structure factor in comparison to the phonon spectrum. 
The spin, dimer, and phonon spectra exhibit a finite size 
gap. It scales with increasing lattice size to a value undis- 
tinguishable from zero, for all our calculations at T = 0. 

Beyond the phase transition, at g = 0.3 (FigQJ), 
the spectrum is different. The renormalized Einstein 
phonons have now joined the Central Peak Branch, form- 
ing a single new phonon branch, which contains the com- 
plete diverging spectral weight. The phonon spectrum is 
affected by the transition only for momenta fairly close 

tO TT. 

The situation is different for larger c^o- Fig. [3] displays 
a density plot with uoo = 1.0 and g = 1 which is close 
to but above the critical coupling. There is a central 
peak branch, which is now separated from the Einstein 
phonons, even above the transition. No softening was 
observed at this larger c^o, even at a large coupling of g = 
1.6. Thus we find a qualitative difference between small 
cjo, where there is a central peak, followed by softening 
beyond g Cl and large c^o, which remains of a central peak 
nature. 

The influence of spins on the phonons is directly visi- 
ble in the dispersion relation of the central peak branch 
at q — TT^ which is is the same as that of the dynamic 
dimer-dimer structure factor 5dd(#, ^) and of the spinon 
branch in the spin spectrum within the accuracy of our 
simulations at all c^o- The dotted lines in Fig.l and Fig. 2 
depict a velocity of v = 7r/2, which is the spin wave 




FIG. 2: Comparison between the dynamical structure fac- 
tor Suu(tt,oo) of dimers SiSi+i and the phonon spectrum 
Sxfau) for cjo = 0.25 and g = 0.23, at N = 256, = 512. 
Whereas the low energy peaks coincide, with a finite size en- 
ergy gap which scales to zero like 1/N, Se>d(tt,^) exhibits no 
structure at the soft phonon peak. The dimer spectrum here 
is almost identical to that of the pure spin model (g = 0). 
The inset shows the spin spectrum (Ss z s z (tt, uo)) as a func- 
tion of k. It reflects the effective increase of the spin coupling 
to about 1.1 J [6( by the phonon interaction, Eq. JT]). 
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FIG. 3: Dynamic structure factor S x (k,Lo) at large bare 
phonon frequency coo = 1 for g = 1, close to but above the 
critical coupling. (N = 256, = 512). 



velocity of the pure Heisenberg chain. This result is ex- 
pected because the operator X{ couples directly to the 
dimer-operator SiSi + i. 

Since the system is critical for g < g c we can also use 
the result of conformal field theory [28[ which predicts 
that the ground state energy Eq(N) can be expressed as 



E (N) =£ (oo) 



6Af 2 



[1 + logarithmic corrections] , 



(4) 

where v is the velocity of the massless modes. Within our 
simulation we find the same velocity v = 7r/2±5% for g < 
g c and coq < 1, at large lattice sizes. When extrapolating 
to the thermodynamic limit, we could not distinguish, 
within our resolution, any logarithmic corrections nor a 
renormalization of the spin wave velocity occuring in the 
large c^o limit [29[ . 
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FIG. 4: Comparison of the phonon spectra S x (tt,cj) at q = tt 
for different system sizes N = 32 (blue) and N = 256 (red) as 
a function of g (cjq = 0.25). While the bigger system exhibits 
a two peak structure with a central peak, the smaller system 
displays softening of the bare phonon mode. 




FIG. 5: S x (q = tt, a;) as a function of spin-phonon coupling 
g for coo = 0.25 and cuo = 1 (inset) at finite temperature 
T = 1/32 (N = 128). 



Systems of small size behave in a qualitatively different 
way. Instead of two branches in S x (k = 7t,cj) they have 
only a one peak structure. The spectrum at N = 32 
in Fig. [H shows a softening of the k = tt phonon mode 
because of the finite size spin gap. At coo = 0.25, large 
systems of TV > 128 are necessary to see the asymptotic 
two peak structure and the new dispersing branch. 

Similar behavior occurs at finite temperatures. Figj5] 
shows density plots of S x (q = 7T,uj) as a function of g 
for uoq — 0.25 and uoq = 1 at T = 1/32. For ljq = 1 the 
situation is similar to the T = limit. For uq = 0.25, 
however, there exits only one single phonon mode which 
softens with increasing coupling. At g ~ 0.3 we find 
a soft phonon mode accompanied by a CPB. Softening 
even without a CPB occurs for c^o = 1 at temperatures 
T>l/8. Approaching the dimerized phase from finite 
temperature, the phonon dynamics ID then resembles 
the 3D case, where for a small ratio of ooo/g a soft phonon 
mode is found, whereas for a bigger ratio the transition 
is governed by a central peak [111 ]. 

Thus we find the following lattice dynamics for the 
T = transition in ID, with a qualitative dependence on 



cjo, which resembles the 3D behavior. As soon as g ^ 0, 
a new narrow phonon branch emerges, centered at zero 
frequency and q = TT. It is clearly separated from the Ein- 
stein phonons at c^o- For big c^o, the Einstein phonons are 
almost not affected by the spin-phonon interaction and 
stay separated from the central peak branch, well beyond 
the phase transition. For small cjo, the Einstein phonons 
soften significantly with increasing coupling before the 
transition takes place. For low enough c^o they soften so 
strongly that they eventually form a single branch with 
the central peak part. Thus for low enough cjo, the T = 
structural phase transition presented here differs from 
its finite temperature counterpart in higher dimension in 
the sense that the phonons soften significantly even in 
the presence of a central peak. 

To conclude, we have directly verified that the quan- 
tum phase transition to the dimerized state is of the cen- 
tral peak type. However, for low uoq <C 1 the Einstein 
phonons soften significantly and eventually join the cen- 
tral peak branch, forming a new phonon branch. 
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